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We investigate matter density perturbations in models of structure formation with or without 
causal/acausal source. Under the fluid approximation in the linear theory, we first derive full 
perturbation equations in flat space with a cosmological constant A. We then use Green-function 
technique to obtain analytic solutions for matter perturbations in a flat A = model. Some incorrect 
solutions in the literature are corrected here. A simple yet accurate extrapolation scheme is then 
proposed to obtain solutions in curved or A / cosmologies. Some general features of these solutions 
are revealed. In particular, we analytically prove that the resulting matter density perturbations 
are independent of the way the causal source was compensated into the background contents of 
the universe when it was first formed. We also use our Green-function solutions to investigate 
the compensation mechanism for perturbations with causal seeds, and yield a mathematically and 
physically explicit form in interpreting it. We found that the compensation scale depends not only 
on the dynamics of the universe, but also on the properties of the seeds near the horizon scale. It 
can be accurately located by employing our Green functions. 



I. INTRODUCTION 



The standard cosmology was lack of a mechanism to produce cosmological perturbations. In order to compensate 
' for this flaw in the standard model, there are currently two main paradigms for structure formation — inflation jl]] 
and topological defects (|] . While the beauty and simplicity of the former appears to have enticed more adherents 
and studies, the latter has proved computationally much more challenging to make robust predictions with which to 
confront observations p|-|T^] . These two paradigms are fundamentally different in the way they generate cosmological 
perturbations. The standard adiabatic inflation produces primordial perturbations on all scales of cosmological interest 
via quantum fluctuations and the causal constraint during inflation, and these perturbations grow over time in an 
uncorrelated manner. As a consequence, the perturbations today can be thought of as simply transfered from the 
initial irregularities that inflation set up, and this transfer function can be easily obtained in the linear theory and 
thus well understood in the literature. On the other hand, topological defects are the byproducts of the spontaneous 
symmetry-breaking phase transition in the early universe, and hence carry energy that was carved out of the originally 
homogeneous background energy of the universe. Therefore due to causality, defects induce perturbations only on 
sub-horizon scales, via gravitational interactions while evolving. This mechanism that prevents the growth of super- 
horizon perturbations is called the 'compensation mechanism'. In addition, due to the certain topology of the defect 
network, the resulting perturbations are correlated and thus non-Gaussian, in contrast to the standard adiabatic 
inflationary perturbations. It then follows that to compute the perturbations in models with defects, we need to 
know the evolution history of defects for the entire dynamic range during which the cosmological perturbations of our 
interest were seeded. This is what makes the computation of defect-induced perturbations so difficult. 

In the literature the power spectra of this kind of models have been investigated using the full Einstein-Boltzmann 
equations. However, the study of the phase information of these perturbations still remains difficult because of 
the limited computation power. Although there have been some detailed treatments for theories with causal seeds 
|T^ . [l5| , we shall in this paper present a simpler formalism, which is an approximation to the full Einstein-Boltzmann 
equations, to provide not only a physically transparent way for understanding the evolution of density perturbations 
in models with source, but also a computationally economical scheme to investigate the phase information of the 
resulting cosmological perturbations. This formalism is parallel to those presented in Ref. [^6| and Ref. |l7| ], but we 
give some modifications to incorporate the inclusion of the cosmological constant and a more detailed treatment for 
the effect of baryon-photon coupling/decoupling. We also note that part of the solutions in Ref. jl6| are incorrect 
due to the incorrect initial conditions and the incorrect assumptions about the form of the subsequent perturbations 
induced by the source (see text later). We shall correct these mistakes and further provide a complete and explicit 
set of analytic solutions for the matter density perturbations. With an accurate extrapolation scheme, these solutions 
become also valid for models with any reasonably chosen background cosmology. The formalism and its solutions to 
be developed here will be completely general and thus suitable for any models with or without causal/acausal source. 
The structure of this paper is as follows. In section O, under the fluid approximation, we first derive in the 
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synchronous gauge the full perturbation equations with source terms, in flat cosmologies with a cosmological constant 
A. This is done by considering the stress-energy conservation of the fluids (II A) and the source (II B), and the linearly 
perturbed Einstein equations (II C). The fluid components considered here are cold dark matter (CDM), baryons (B), 
and photons (7), and we employ the baryon-photon tight-coupling approximation to derive the perturbation equations 
before the last-scattering epoch. In this context, we also investigate the role of the so-called stress-energy pseudotensor 
(II D). The initial conditions of these perturbation equations are discussed (HE), and we use the approximation of 
instantaneous decoupling to deal with the decoupling of photons and baryons at the epoch of last scattering ( |II F| ) . 
We then numerically justify the accuracy of this formalism in the context of standard CDM models, by comparing 
its results with those of the full Einstcin-Boltzmann solver |13| (II G). Within reasonable ranges of cosmological 
parameters, our approach provides satisfactory precision at greatly reduced numerical cost. 

In section III, we derive the matter perturbation solutions of the equations presented in section [n| The perturbations 
of radiation and matter are first divided into two parts: the initial and the subsequent perturbation s. W ith some 
change of variables, these equations are then ready to be solved by the Green- function technique ( III A ). With 
this technique, we find the exact solutions on scales much larg er or much smaller than the horizon size, namely 
the super-horizon or the sub- horizon solutions respectively ( III B ). Some degeneracy among the G reen functions for 
the matter perturbation solutions is then found and used to reduce their effective number (QIC). With this great 
simplification, solutions on intermediate scales are then easily obtained by an accurate interpolation scheme based 
on the well-known standard CDM transfer function ( |III D| ) . We also discuss the effect of baryons ( III E ). A simpl e 
and accurate extrapolation scheme is then introduced to obtain solutions in the K ^ or A ^ cosmologies ( IIIF ), 
where K is the curvature of the universe (see Appendix |X|) . All our Green- function solutions are numerically verified 
to high accuracy. 

In section IV, we use our Green- function solutions to investigate some important properties of cosmological matter 
density perturbations. We first demonstrate the relation between our solutions and the standard CDM transfer 
function (IV A). We also prove that in models with causal source, the resulting matter perturbations today are 
independent of the way the source energy is initially compensated into the background contents of the universe 
( IV B ) . Finally we use our Green- function solutions to study the compensation mechanism and the scale on which it 
operates (IV C). We find that this compensation scale is determined not only by the dynamics of the universe, but also 
by the properties of the source near the horizon scale. Once the detailed features of the source near the horizon scale 
are known, this compensation scale can be accurately located using our Green functions. A summary and conclusion 
is given in section [v|. In appendix |A|, we define the convention of some notations used in this paper, and present for 
reference the solutions for the dynamics of various background cosmologies, including the consideration of non-zero 
curvature and a cosmological constant. 



II. SYNCHRONOUS GAUGE PERTURBATION THEORY 



In this section, we derive the linear evolution equations for cosmological perturbations. To calculate the density 
and metric perturbations, we model the contents of the universe as perfect fluids: radiation (photons and neutrinos) 
and pressureless matter (CDM and baryons). We shall use the photon-baryon tight-coupling approximation until 
the epoch of last scattering, at which we assume instantaneous decoupling, also taking into account the effect of 
Silk damping due to the photon diffusion. After the decoupling, the baryonic perturbations originating from the 
perturbations of the photon-baryon coupled fluid are then merged linearly into the CDM content. In scenarios with 
causal seeds, the radiation and matter fields are assumed to be initially uniform, and then perturbed by the causal 
seeds after they are formed. The radiation, matter, and causal seeds are assumed to interact only through gravity, 
meaning that their stress-energy tensors are separately covariantly conserved. 

We shall work in the synchronous gauge, in which the perturbations h^ v to the spacetime metric obey the 
constraint ho^ — 0. Throughout this paper, we use a signature ( — h H — h) for the spacetime metric, and units in which 
% = c = kg = 1. Thus the perturbed flat Friedmann- Robertson- Walker (FRW) metric is given by 

ffoo = -a 2 (rj), 9tj = a 2 {r])[S tj + ^(77, x)]. (2.1) 

We shall work in the linear theory, requiring \hij\ <C I- Greek alphabet will denote the spacetime indices (e.g. 
fx = 0, 1, 2, 3), and mid-alphabet Latin letters the spatial indices (e.g. i = 1,2,3). Although the synchronous gauge is 
sometimes criticised in the literature due to its residual gauge freedom, it is still well suited to models in which the 
universe evolves from being perfectly homogeneous and isotropic. In such models, all perturbation variables can be 
initially set to zero (before the causal seeds are generated), and this is normally referred to as the 'initially unperturbed 
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synchronous gauge' (IUSG) [[16]]. It possesses no residual gauge freedom. Thus the Einstein equations are completely 
causal in IUSG, with the values of all perturbation variables at a given spacetime point being completely determined 
by initial conditions within the past light cone of the point. One example of such models is the cosmic defect models, 
which have been of most interest in the study of models with causal seeds. 

we 



In section II A, we derive in the IUSG the conservation equations of radiation and matter fields. In section [IB 



consider the conservation of source stress energy. In section II C, we derive the linearly perturbed Einstein equations. 
Then, in section II D, we employ the concept of stress-energy pseudotensor to investigate the internal energy transfer 
among various fields. In section II F, we describe the approximation of instantaneous decoupling of photons and 
baryons at the epoch of last scattering. In section [I G , we numerically verify the accuracy of our formalism for the 
standard CDM model, in comparison with the results from CMBFAST [n~3j , a fast Einstein-Boltzmann solver. 



A. Stress-energy conservation of radiation and matter fields 

The contents of the universe are considered as perfect fluids, whose energy-momentum tensors have the form 

T nu = (Pn + Pn)Kunu + PnKi with K u nh = -1- (2.2) 

Here p N , p N , and u N are the density, pressure, and four- velocity of the Nth fluid respectively. In the homogeneous 
background, we have = (a -1 , 0), which implies that Su° N = to first order. We thus define the velocity perturbation 
as v l N = a8u l N , i.e., 8u N = (0, ~v N /a). The equation of state and the sound speed are defined respectively as 

n N = P -^, cl= 6 ^. (2.3) 

Pn 0p N 

Consequently, the covariant conservation of stress energy for each fluid Tj}" v = gives jl6| 

S N + (1 + k w )(V • v N + -h) + 3-(4 - k n )S n = 0, (2.4) 
2 a 

a c 2 
v N + -(l-3c%)v N + —^V6 N = 0, (2.5) 
a 1 + k n 

v^ + -(l-3c 2 „)v^=0, (2.6) 
a 

where S N = Sp N /p N , h = ha is the spatial trace of h^, and we have decomposed the velocities as v„ = v N + v^, 
with V x v|i = and V • = 0. 

In the regime of photon-baryon tight coupling, we have only two main fluids: the CDM component and the tightly- 
coupled photon-baryon fluid. They will be denoted as N — c, 7B respectively, and discussed separately as follows. 
Note that we have ignored the neutrinos in the radiation. 



1. CDM fluid 



We first consider the CDM fluid, i.e. N = c. With n c 



conservation (2.4)-( |2.6| ) become 



for pressureless matter, the equations of stress-energy 



V- v c 



-h, 



0. 



(2.7) 



As we can see, any perturbations in the CDM velocity will decay as a -1 . Thus we can simply choose v c = in the 
IUSG. Once v c = 0, it will remain so as there is no linear gravitational source. As a consequence, the CDM obeys a 
single nontrivial conservation law resulting from equation (2.7) 



h + 25 c = =^> h = -26 c , 
where the second equation results from the initial condition h = S c = 0, as required by the IUSG. 



(2. 
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2. Photon-baryon tightly coupled fluid and its photon component 



For the tightly-coupled photon-baryon (7B) fluid, we have 

v 7 b = v 7 =v b , p 7 B=P 7 , P 7 b = P 7 +Pb- 

Thus we can define 

R = ^Pb _ 3pB 
5p 7 4p 7 ' 

where the second result comes from the fact that p 7 cx a -4 and pb oc a -3 . Definitions (2.3) then give 

1 , 1 

K 7 B 



(2.9) 



(2.10) 



(2.11) 



3 + 4i?' 78 3(1 + R)' 

With these results, the equations of stress-energy conservation for the 7B fluid can be obtained from equations 
(O-O): 



4 + AR 
3 + 4i? 



(V • v 



7 b 



i . a 
6c) + ~ 



R 



a (I + R){3 + 4R) 



<5 7 b = 0, 



V 7 B + 



a R 



al + R 



v 7 b 



3 + 4i? 
12(1 + i?) 2 
1 d R 



V(5 7 b = 0, 



V ^ + aTTR V ^ = °- 



(2.12) 
(2.13) 
(2.14) 



In cosmological applications, such as CMB anisotropics, we are more in tere sted i n the photon perturbations rather 
than the perturbations in the jB fluid. Therefore by using equations (^j]) and ( 2.1C ), we can extract the photon 
component from the above equations to yield [fL8| 



4 c 

-V • v r 5 C 

3 3 c 



0, 



a R 



1 



v r + - 



al + R 



V r 



4R 



ws r = 0, 



(2.15) 
(2.16) 



where we have ignored neutrinos in the radiation so as to replace the subscript 7 with r. The velocity can then be 
eliminated to yield a single second-order equation: 



Sr 



—^(4- -s c ) 
i + ir 3 ' 



1 



3 + 3i? 



V 2 (5 r = 0. 



(2.17) 



We note that although the photon velocities are missing in this equation, they can be recovered at any given moment 
using equation ( p. 15 ). 

An alternative presentation of equations ( 2.15| ) and (2.17) is via the entropy perturbation s. It is defined as the 

(2.18) 



fluctuation in the number of photons per dark matter particle 

3 



s = -S r ~ S c 
4 



Thus equations ( p. 15 ) and (2.17) can be rewritten as 

s = -V • v r , 
R 



1 



1 + R 3 + 3i? 



V 2 (s + 5 C ) 



(2.19) 
(2.20) 



As we shall see, 5 T can only have a white noise power spectrum on super-horizon scales. From equation (2.16), this 
implies a k 2 power spectrum in v r on these scales. Adding the fact that the entropy fluctuation s starts from zero on 
super-horizon scales due to the fixed number of dark matter particles per photon, it then follows from equation (2. IE) 
that both s and s have a fc 4 fall off outside the horizon. Therefore in numerical simulations, as long as the initial 
horizon size is smaller than the scales of our interest, we can simply set s = s = as part of the initial condition. 
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B. Stress-energy conservation of the source 

The causal source we shall consider is weak, so it will appear only as first-order terms in the perturbed Einstein 
equations. Thus in the linear theory we are considering here, they can be treated as being stiff, meaning that their 
evolution depends only on their own self-interactions and the background dynamics of the universe, but not on their 
self-gravity or on the weak gravitational field of the inhomogeneities they produce. This assumption will enable us to 
separate the calculation of their dynamics from that of the inhomogeneities they induce, allowing us to evolve them 
as if they are in a completely homogeneous background. Since the source is stiff, its energy-momentum tensor 
need only be locally covariantly conserved with respect to the background: 

e 00 ,o + -e+ = e 0M , (2.21) 

a 

6 oi ,o + 2% 0l = <->,,, , (2.22) 

where 0+ = 600 + @u- 

Another important aspect of cosmic structure formation with causal seeds like cosmic defects is the fact that the 
sources, formed at very early times, will ultimately create under-densities in the initially homogeneous background, 
out of which they are carved. This is a direct result of energy conservation in the universe, and is normally termed 
'compensation'. We shall discuss this issue in more detail later. 



C. Linearly Perturbed Einstein equations 

At first we have ten Einstein equations 



or equivalcntly, 



R„ v = 8nG(T^ - l -g^T) + Ag„ v , (2.23) 

G>„ = - \g^Rs = SttGT^ - Ag^, (2.24) 

where is the Ricci tensor, G is the gravitational constant, T = g^T^, A is the cosmological constant, G^ is the 
Einstein tensor, and i?s = g^ v R^u is the scalar curvature. Linearly perturbing the above equations, we obtain 

SR^ = 8nG(5T^ - h^ V rs T rs - l V ^ V rs dT rs + ^/^V 9 ^) + a 2 Ah^, (2.25) 



or equivalcntly 



SG^ = 8irG6T^ - a 2 Ah^, (2.26) 



where 

ST^ = 0^ + a 2 Y,(h^ r T r N „ + ^s5T s Nu ). (2.27) 

N 

A closed set of the ten linearly perturbed Einstein equations are then 



25Ro = -h - % = +3 (^) ^2(l + 3c 2 N )n N 6 N + 8nGG+ , (2.28) 



N 



2 < SRij 



2 

a la' 



+ 

a \a 



2SGoo = hij,ij -\v 2 h + 2-h = &[-\ ^ £l N 6 N + 167rG6 o , (2.30) 



hij \ = hij + 2-hij - \7 2 hij - \h^j + ^S i3 \7 2 h 

+hik,kj + hjkM — -^5ijhki,ki — lQnGQij + a Ahij , (2.29) 



3 



26G M = hij d - h 4 = -6 ( - ) ^ C 1 + + 16^G9 0J , (2.31) 







•a 


) 








a 








a 



N 
2 



N 
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where the traceless parts are defined by R^ = Rij — 5ijR k k/3, and similarly for hij and Oy ■. The prime over the sum 
in equation (2.31) indicates the sum over all fluids except CDM. We note from the above results that in the IUSG 
the cosmological constant A does not appear as extra terms in the perturbation equations except in (2.29), the Hj' 
component. 

Within the photon-baryon tight-coupling regime, the above perturbation equations simplify as: 



h- a -h = 



-3 



-) 1(2 



R)Q r S r + n c <y + 87TG6+ 



2-L 



1 1 2 ~ 2 - 



a • 
2-h = 6 

a 



[n c S c + (1 + fl)ft r( J r ] + 167rG6o 



H/^ = -8 I - I (1 + E)n r «J + 167TG6 



CH 



(2.32) 
(2.33) 
(2.34) 

(2.35) 



We note that if the source obeys the covariant conservation equations (2.21) and (2.22), then equations ( [2.34 ) and 
( |2.35 ) are preserved by equations (2.32). 

In the standard CDM model where the source is absent, equation (2.32) can be greatly simplified on super-horizon 
scales {krj <§C 1) in the radiation or matter era: 



5c + -5, 



1- 2(2 + R) 



V 2 

v 



6 C = 0, in radiation era, 



;S C = 0, in matter era. 



(2.36) 
(2.37) 



Since R = 3f2Bo a /4^ c o a cq by definition, we know R <C 1 deep in the radiation era. Thus the above equations both have 
a growing mode <5 C oc rj 2 . This result has an important implication for numerical simulations of structure formation 
with causal sources. In this case, if numerical errors appear as white noise on super-horizon modes k < 1 / 77, then 
they will have a growing behavior S(k) = 47rfc 3 'P(&;) oc k?rf. For the horizon crossing mode k ~ l/?7, this becomes 
S(k) cx T] jL7|. This means that although energy conservation together with causality should forbid the growth of 
perturbations on super-horizon scales, any numerical errors seeded from early times would induce a spurious growing 
mode on these scales. To overcome this problem, one needs to perfectly compensate the source energy in the initially 
homogeneous background. In the following section, we shall discuss one of the methods that can achieve this. 



D. Stress-energy conservation of the pseudotensor 

The concept of the stress-energy pseudotensor in an expanding universe was first remarked in this context by 
Veeraraghavan and Stebbins [^6[ , and further investigated by Pen, Spergel and Turok fl7j| . To introduce this concept, 
we start from a perturbed Minkowski space = rj^ + h^, where the Bianchi identity V^G^ = leads to an 
ordinary conservation law c^G^i) — at linear order in h^ u . Adding the fact that the Einstein equations give 

= SirGT^u — G^ v i n \) where G^ni) is the sum of non- linear terms in /i M „, we see that the right-hand side of 
this equation provides an ordinarily conserved tensor, the stress-energy pseudotensor. 

The generalization of this result to an FRW model is straightforward, with only the corrections due to the expansion 
of the universe. Moving all these corrections (derivatives of the scale factor) to the right-hand side of the Einstein 
equations while keeping only the linear terms in h^ u , we obtain a pseudo-stress-energy tensor t^ v = G i1u <i)/S>'kG: 

T o° = 7T7;(-) [n c S c + (l+R)n i 5 I }--^—-h + e 00 , (2.38) 
87rG \a J oirG a 

1 I a 



2ttG 



T0i = ~^ - (l + iOfirvJ + eoi, (2.39) 



1 ( a\ 2 „ „ 1 a,~ 2, 



T - = 5 -^G UJ QA ~ ^Ga {H - ~ 3^ + 9 « ' (2 - 40) 
This tensor obeys an ordinary conservation law v = according to the Einstein equations, or equivalently 



G 



no 



i0 



Jli 



(2.41) 
(2.42) 



This is not a fundamentally new conservation law, but it describes the interchange of energy and momentum among 
the different components in the universe, i.e. the radiation, matter, and the source in our case. This description 
appears to be physically more transparent than the original Einstein equations. 

Another advantage of invoking this formalism is that it is easier for numerical simulations to specify the initial 
conditions and to maintain proper compensation on super-horizon scales. As we shall explain later, ry can only have 
a white-noise power spectrum on super-horizon scales. Thus integrating equations (2.41) and (2.42) over time shows 
that Too has a fc 4 power spectrum and that Tbi has a fc 2 power spectrum. Therefore, as long as the horizon size at 
the beginning of the simulation is smaller than the scales of our interest, we can set Too = T oi — as the initial 
condition, allowing for perturbations to grow only inside the horizon and for Too to fall off as fc 4 outside the horizon. 
For simulations of structure formation with causal source, a check of Too °t k 4 on super-horizon modes will tell us 
whether or not the compensation is well obeyed. 

To make use of the pseudo-stress-energy tensor formalism in the study of cosmological perturbations, we combine 
the conservation equation for radiation ( 2.20| ), the definition of pseudoenergy (2.38), and one of the alternative Einstein 
equations using the pseudo-stress-energy tensor (2.41), to yield a convenient closed set of equations: 



R 



1 



1 + R 3 + 3R 



V 2 (s + 4) 



S c = 4ttG-(t o - Goo) - - 
a a 



-n c + 2(i + R)n r 



6 C + 2(1 + R)n r s 



TOO = Go 



1 fa 
2^G [a 



(1 + R)Q T s 



(2.43) 
(2.44) 

(2.45) 



Here we have used equations (2 
the results in Ref. 
formalism. 



(2. IS), (2.19) and (2.39) to eliminate h, S r 



and Toi respectively. By analogy to 
here we have built both the pseudoenergy tqq and the entropy fluctuation s into the above 



E. Initial conditions of causal models 



As required by the IUSG, all perturbation variables are zero before any mechanism of structure formation starts to 
act on the initially homogeneous and isotropic universe. In causal models, causality also requires that local physical 
processes can never induce correlated perturbations on scales much larger than the horizon. Therefore, when the initial 
irregularities of the universe are first formed (e.g. via the formation of cosmic defects, or the presence of inflation), 
the spatial part of 0^, t^ v and h^ v can only have white- noise power spectra on super- horizon modes — their spatial 
perturbations being uncorrelated on scales larger than the horizon size fl6[| . The same applies to <5/v and therefore h. 
It then follows from equations ( 2. 16 ), ( 2.22 ) a nd (2.42 ) re spec tively that the power spectra of v r , Goi and Toi all fall 
off as k 2 outside the horizon. From equations (2.19) and (2.41), we also have the spectra of s, s and Too proportional 
to fc on these scales, as previously discussed. As a summary, we have for super-horizon modes k -^l/f] that 



$c ; <^r ; h, flij , Qij , Tij oc k , 

v r ,Goi,Toi oc k 1 , 
s,s,t 00 oc k 2 , 



(2.46) 
(2.47) 
(2.48) 



where 'oc fc™' means the power spectrum is proportional to k 2n . 

Since the production time of the initial irregularities is normally so early that the horizon size rji at that time is 
much smaller than the cosmological scales k~^ s of our interest (i.e. fc C os??i *C 1), the above conditions can be regarded 
as general initial conditions for all scales of cosmological interest. If we require fc CO s??i 1 in our analysis, we can 
simply choose 



Vl = Qoi = T i = S = S = Too = 0, 

as the initial conditions, because their power spectra all decay as either fc 2 or fc 4 outside the horizon. 



(2.49) 
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With such a choice, we can see from equation ( 2.38| ) that there is still freedom for the choice of 8 C , S T and h into 
which to compensate 0oo when 0oo was first formed. Nevertheless, as we shall analytically prove later, no matter 
how 600 was compensated into the background contents of the universe when the causal source was first formed, the 
resulting matter density perturbations today would be the same. We note that this was first numerically observed 
in Ref. ]l7j j, and here we shall provide a thorough interpretation to it using our analytical solutions to be obtained 
later. We also note that none of the above arguments will hold if the initial perturbations are seeded in an acausal 
way, which is nevertheless not of our current interest. 



F. approximation of instantaneous decoupling 

One thing we have not included in our formalism is the treatment at and after the decoupling epoch r/d- Before this 
epoch, photons and baryons are assumed to be tightly coupled, forming a single 7B fluid. At the decoupling epoch 
?7d, baryons and photons are assumed to be instantaneously decoupled from each other, so that <5 7 and 5b evolve 
separately afterwards. A numerical fit to the redshift of the decoupling epoch is jfj) 

h = 0.313(n m0 h 2 )- 0A19 [1 + 0.607(a„o/i 2 ) (574 ] , (2.51) 
b 2 = 0.238(fi m0 /i 2 ) - 223 . (2.52) 

Although this is the result for the decoupling epoch of baryons and there is another fit for that of photons, these two 
epochs — the recombination of baryons and the last scattering of photons — coincide approximately in the absence of 
subsequent reionization pp| , pl[ . 

In addition, the photons and baryons are not in fact perfectly coupled, and this leads to the diffusion damping of 
photons and Silk damping of baryons |p2|| during the decoupling epoch. To model these effects, we apply damping 
envelopes to both (5 7 and <5b at the decoupling epoch Zd, i.e. 

S N (d) = S N{d) V N (k), -/V = 7, B, (2.53) 

where the tilde indicates the Fourier transform of a quantity and k is the wave number. The photon diffusion damping 
envelope can be approximated by the form 

£> 7 (fc) ~ e - {k/k - )m \ (2.54) 

where 



/c 7 I 2 



Mpc -1 I 7T 



arctan 



7T /F 2 ^ P2/P1 



2 V f\ 



(n B0 h 2 r 2 



pl/p2 



Fi , (2.55) 



m 7 = 1.46(a n0 /i 2 ) 0303 (1 + 0.128 arctan {in [(32.8ft B o/i 2 r°- 643 ] }) , (2.56) 

pi = 0.29, (2.57) 

p 2 = 2.38(a n0 /i 2 ) - 184 , (2.58) 

Fx = 0.293(ft m0 /i 2 ) - 545 [1 + (25.m m0 /i 2 r a648 ] , (2.59) 

F 2 = 0.524(fi m0 /i 2 ) - 505 [1 + (10.5fi m0 /i 2 )- - 564 ] . (2.60) 

Silk damping for the baryons can likewise be approximated as pl| 

V B (k) ~e- (fe/fcs)ms , (2.61) 



where 



fc S _ 100/0 ^0.398,o ,2,0.487 l + (96.2anO^- !U S 



Mpc~^ ^ ^^^^^ T^ag&pacr ' (2 ' 62) 

, n (»B0fe 2 )-° 0297 ffl m 0fe 2 ) 0282 

WS = L4 ° i + (78inBofr a )-°-«"« ■ (2 ' 63) 
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In some scenarios with causal sources, the damping envelopes (2.54) and (2.61) may depart from the form of 
exponential fall-off here to a power-law decay towards smaller scales. This is due to the survival of perturbations 
which are actively seeded during the decoupling process. For example, in models with cosmic strings, the departure 
appears on scales smaller than of order a few arc- minutes (i.e. the multipole index I > 3000) [^3|. Certainly this 
is beyond the scale range of our interest. Moreover, since the decoupling process is relatively a short instant in the 
entire evolution history of the perturbations, the contribution from these survived small-scale perturbations should 
be relatively small. Adding the fact that we expect the post-decoupling contribution in the perturbations seeded by 
defects to have a power-law fall-off on small scales due to a certain topology of the source || , the small-scale power 
in the final perturbations is likely to be dominated by this post-decoupling contribution, rather than the primary 
perturbations (those seeded before and during the decoupling, whose power spectrum first exponentially decays and 
then turns to a power-law fall-off). Therefore, on the scales of our interest, the damping approximation employed 
here should be still appropriate for models with cosmic defects. 

No w we consider the evolution of <5 7 and Sb after the decoupling epoch z&. From the energy conservation law 



fl2.4|) -(2.6), we have for the baryon perturbations 



a 



S b -S c + -(Sb-5 c )=0. (2.64) 



a 



This implies (Sb — S c ) cx a -1 , meaning that the evolution of Sb and <5 C will soon converge to the same behavior. 
We also know that matter perturbations grow as rf 1 in the matter era so that [<Wd) — ^c(d)] is relatively small when 
compared with either Sbo or S c q. As a consequence, in the calculation of <5bo and S c q to linear order, it is appropriate 
to combine Sb and S c at the decoupling epoch Zd as 



^B0<5B(d) + ^c0<5c(d) _ 3f2BO<5 7 (d)f b/4 + ^c0i5c(d) 



"m(d) = Fi r~F> = o ~ ' (2.65) 



and the same for their time derivatives. Then we have only two fluids after the decoupling: the photon fluid (f2 7 ) and 
the matter fluid, which is linearly combined from the CDM and baryon fluids (f2 m = Q c + Qb)- Eventually we can 
take the matter perturbations at the present epoch to be S c0 = Sbo = S m o- 



To sum up, we first evolve the CDM and 7B perturbations up to the decoupling epoch Zd given by ( 2.5C ), noting 



that our formalism extracts the photon component 7 from the 7B fluid. We then apply damping envelopes to <5 7 (d) 



and <5 c (d), as illustrated by equation (2.53), to account for the photon diffusion and Silk damping. <5 m (d) is then 



obtained by linearly combining <5 c (d) and <Wd); as shown in equation ( 2.65 ). Finally we carry on the evolution of S m 
and 5 r from the epoch Zd to the present, using our previous perturbation equations with R — and the subscript 'c' 
replaced by 'm'. 



G. Accuracy for the standard CDM models 



To verify our scheme for evolving cosmological perturbations, we first calculate the CDM transfer function in the 
context of the adiabatic inflationary CDM model: 



T c (k,T) ) 



Sc(k, m )S c (0,0) 
S c (k,0)S c (0,r) y 



(2.66) 



where 770 is the present conformal time. To this end, we employ equations (2.43), (2.44) and (2.45) in the absence 
of the source terms, and the approximation of instantaneous decoupling described above. We start the evolution in 
the deep radiation era when O m -C O r « 1, R -C 1, and 771 <C 1/k for a given mode k. In this case, one choice of the 
initial conditions is 



0, 5 C 



V?, 



too 



1 

ttG' 



(2.67) 



Figure [l] shows our results for the CDM transfer functions T c (k,rj ) at the present epoch in different cosmologies, 
together with the results obtained from CMBFAST |13|. It is clear that they agree very well. The discrepancy of the 
two reaches its maximum of about 5% at the scale k ~ 1/iMpc -1 in the open model with fi c o = 0.15 and £1bo = 0.05. 
We have also checked our results against those in Ref. p4[, and they are in agreement again within a 5% error. In 
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k(h Mpc" 1 ) 

FIG. 1. Comparison of our CDM transfer functions at the present epoch T c (k,rjo) with results obtained from CMBFAST 
[Jl3|. On the top are results in flat models with a cosmological constant (i.e. Qao + ^co + ^bo = X). At the bottom are results 
in open models without a cosmological constant. Results using our formalism are plotted as solid lines, while the results from 
CMBFAST are plotted as dashed lines. We have used h — 0.7 throughout. The mass fraction of Helium-4 Yh c = 0.24 and the 
number of neutrino species N v = 3.04 have been used in obtaining the results from CMBFAST. 
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addition, from the bottom curves in Figure H, we notice the oscillations resulting from the photon-baryon coupling 
before rj^ in cosmologies with high baryon fractions fiso/^mo- 

Next, we calculate the radiation transfer function at the decoupling epoch, since the radiation perturbations at this 
epoch will appear as the intrinsic CMB anisotropies. We define this transfer function as: 

T v (k,r] d ) = - , (2.68) 

S c (k, 0)^(0,770) 

where we normalize the radiation perturbations at r/ d to both the amplitude of the super-horizon CDM perturbations 



today and the initial CDM power spectrum, as we did for T c (fc,?7o) (sec cq. [2.66]). This definition will enable us to 
verify not only the scale dependence of the evolution of perturbations, but also their normalizations. Figure [2] shows 
our results, again as a comparison with the results from CMBFAST. We see that although the scale dependence of 
our results is slightly different from that of the CMBFAST results, the overall normalisation appears to be quite 
accurate. The sideway shift of the oscillatory peaks in our results when compared with the peaks from CMBFAST 
has a maximum of about 5% in the flat model with Q c q = 0.95 and f2eo = 0.05. This discrepancy results naturally 
from the instantaneous-decoupling approximation in our formalism. As a result, despite the small inaccuracy, our 
formalism provides a much more numerically efficient way than the full Einstcin-Boltzmann scheme in calculating the 
density perturbations. 



III. SOLUTIONS OF MATTER PERTURBATIONS 



A. Decomposition of perturbations 



We first consider density perturbations about a flat FRW model with a cosmological constant A, which are causally 
sourced by an evolving source field with the energy- momentum tensor MJ/ (x, 77) . As seen in the previous section, 
with the photon-baryon tight coupling approximation in the synchronous gauge, the linear evolution equations of the 
radiation and CDM perturbations can be given by equations (2.43), (2.44) and (2.45), which are derived from equations 
(2.20), ( 2.38| ) and (2.41). This set of equations has the advantage in controling the initial condition for numerical 
simulations, as well as understanding the law of stress-energy conservation. For analytic simplicity, however, we shall 
drop the use of Too m this section, and employ equations (2.17) and (2.32) to form an alternative set of evolution 
equations for density perturbations: 



8\ 

3 

: a ■ 



R 



1 

3/dv- 
2 \a 



3^ c) 3(1 + R) 



V 2 <5 r = 0, 



(J^J [n c 5 c + (2 + R)n T S T ] = 47rG6 H 



(3.1) 
(3.2) 



We note again that the cosmological constant A affects only the background dynamics (i.e., the evolution of the scale 
factor a), but does not contribute extra terms in the above perturbation equations. After the decoupling epoch rj^, 
the treatment is essentially the same as that introduced in section [IF. We have numerical verified in the context of 



the adiabatic inflationary CDM model that the set of equations (3.1) and (3.2) and the set of equations (2.43), (2.44) 
and (2.45) indeed give identical transfer functions of density perturbations, with a numerical discrepancy of less than 
0.1%. 

Assuming that the causal source was formed at some initial time 77; and then evolved to the current time 77, it proves 
useful to split the source-seeded linear perturbations into initial (I) and subsequent (S) parts E3|: 



S N (x, 77) = S l N (x, 77) + Sff(x, 77) , N = c, r. 



(3.3) 



The initial perturbations <5]y(x, 77) originate from the source configuration at 77;, while the subsequent perturbations 
<5)y(x, 77) are actively and cumulatively seeded by the later evolution of the source at each 77, where 771 < 77 < 77. This 
is equivalent to having the initial conditions 



8%(m) = s%(m) = 0. 



(3.4) 
(3.5) 



Because the source induces isocurvature perturbations, 5 : (x, 77) must compensate S s (x, 77) on comoving scales |x— x'| > 
77 to prevent acausal perturbation growth on super-horizon scales. One of the aims of this paper is to show analytically 
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-0.5 1 



10" 



10 



10" 



k(h Mpc -1 ) 

FIG. 2. Comparison of our radiation transfer functions at the decoupling epoch Tr(fc,7?d) with results obtained from CMB- 
FAST jl3| . On the top are results in flat models with a cosmological constant (i.e. D.ao + ^co + ^bo = 1)- At the bottom are 
results in open models without a cosmological constant. Our results are plotted as solid, dot-dashed, and dotted lines, while 
the CMBFAST results are plotted as dashed lines. Note that this transfer function has been n orma lised to both the amplitude 
of the super-horizon CDM perturbations today and the initial CDM power spectrum (see eq. 2.68]). 
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how this compensation mechanism can be achieved. Now we can solve the system of equations (3.1) and ( p.2[ ) by 
employing the integral equation with Green functions: 

(3.6) 
(3.7) 



5 1 N (^V)=J2 d 3 x'g? N '(X;r ll r ]i )5 N ,(x' im )+ U 3 x> G? N ' (X; r,, m )6 N ,(x! , m ), 
N , J J 

S%(x,r)) = 4ttG J" dfj J dV^fliMje+fx',,)), 



where X = |x — x'|. The easiest method of obtaining the Green- function solutions is to go to Fourier space and solve 
the resulting homogeneous system of ordinary differential equations with appropriate initial conditions. Since the 
Green functions depend only on the modulus of X = |x — x'|, it follows that their Fourier amplitudes must depend 
only on the modulus of k. Thus we have 



? N (k, t}) = [Si N \k; V , m)8 N , (k, + G? N '(k;r ] ,r ]i )5 N ,(k,r h ) 

N' 



6 s N (k,ri) = 4wG / £ Ws (fc; 77, 7))9+(k, 77)^77. 

Jr/i 



(3.8) 
(3.9) 



We notice that equation (3.9) is different from the form in Ref. JT(|, where the authors identified our Q Ns as Q. 



Nc 



This identification is incorrect, because Q and G% c have different initial conditions, as we shall see. 

For simplicity, we assume no baryons and therefore set R = for now, and shall relax this constraint later. With 
t he chang e of variable y = 1 + Arj/2 wher e A = 2{\ [2 — l)/r] eq (leading to a/a cq = y 1 — 1), and with the formalism 
( |3.8|) and (3.9), we can rewrite equations (3T) and (|3.2| ) in Fourier space as 



~ 4 ~ 4fc 2 ~ 
y 3 3A 2 



(l - y z )9 c " - 2 y g c 



i2gyg c 
i-y 2 



(3.10) 
(3.11) 



where a prime represents a derivative with respect to y, Q c = Q1 



cN 



QZ N or g c 



and Q r 



to equations (3.S) and (3.9), the initial conditions (3.4) and (3.5) now become: 



G{ , G^ N or Q rs ■ According 



net 



= Ql r = Qf = 1 at r? 

~ 3 ~ 

g cs = -g rs = i at 77 

4 



fj, 



(3.12) 
(3.13) 



with all the other Green functions and their time derivatives vanishing. There are three things we should notice here. 
First, it is required that (k; rj, r/i) = for r\ < 77;, and that G Ns (k; 77, 77) = for 77 < 77. Second, the Green functions 
Gf N ' only describe the time dependence of the homogeneous version of equations (3.1) and (3.2), while the Green 
functions G Ns are, by the conventional definition, the true Green functions used to solve the inhomogeneous equations 



(3.1) and (3^2). Finally, since there are only four variables in equations ( p.l0| ) and ( 3.11 ) (i.e. G°, G r , Q and Q ), 
there must exist some dependence among the five sets of Gr een fu nctio ns (i.e. G\ c , Gi, G^ C : G2 1 an d G Ns )- This 
dependence can be observed from the initial conditions ( 3.1 2| ) and (3.12), which yield 



G Ns = Gi 



o^2 • 



(3.14) 



In Ref. [^6|, the authors ignored the fact that G ls = 4/3 in the initial condition ( |3.13j ). This ignorance led to the 



absence of the second term in equation (3.14) (and thus the identification of G Ns — G2°)i an d consequently the 



incorrect solutions of Green functions in their final results. Based on equations (3.10) and (3.11) with the initial 
conditions ( 3.12| ) and ( [3. IS ) , in the following subsections we shall analytically derive a complete set of Green- function 
solutions for the matter perturbations, which will then be numerically verified. 
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B. Super- horizon and sub- horizon modes 



Under the limit krj ^§> 1 or krj -C 1, the ratio Q T jQ c will approach a constant (see below), so that equation ( 3.11 ) 
becomes the associated Legendre equation, with solutions composed of the associated Legendre functions P% (y) and 



Qa (y), where fi = y 12Q T /Q C . We shall use subscripts oo and to denote solutions in the limits krj ^> 1 and krj <C 1 
respectively. For simplicity, we shall denote both fj and r/i as fj in the following solutions. 

1. kr/ 3> 1: When the wavelengths are much smaller than the horizon size, the radiation oscillates many times per 



expansion time and its effect is therefore negligible. By setting Q T /Q C = 0, equation (3.11) can be solved as 

J&fo, fj) = E(fj)P°(y) + F(f,)Q° 2 (y), (3.15) 

where E(fj) and F{fj) are functions of fj. This gives the sub-horizon solutions. 
2. kr/ <C 1: When the wavelengths are much longer than the horizon size, we have Q T /Q c — 4/3 as the consequence 



of zero entropy (see eqs. | 2.18 | and | 2.49 |), giving \i — 4. Thus equations ( 3.10 ) and ( 3.11 ) yield 

4 „ 

G T (v, V) = ^Go( r h V) + a i(v ~ V) + A ) 



G c ( V ,fi) = G(r))P 2 - 4 (y) + H(fj)Qt(y) 
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» Qt(a;)P 2 - 4 (y) - P^^Q^) A# + 2 ai (x - Vi ) 
, Q^i^x) - P 2 - 4 (.x)Qf (x) A(£ 2 - l) 2 



(3.16) 



(3.17) 



where a; and /3j are constants, and G(f)) and H(fj) are functions of fj, all determined by the initial conditions. 
These are the super-horizon solutions. 



Combined with the initial conditions ( 3.1 2| ) and (3.12), equations (3.15) and (3.17) can be solved to yield the following 
results. For clarity, we shall denote y = 1 + Afj/2 in Q Ns and yi = 1 + Ar/i/2 in Q^ N both as w : 



Woo 

/-">CC 

yioo 

p>cc 

/ice 
^20 

yioo 

P> cr 

yio 

t/2oo 
^20 



— (™ 2 -l)<j(3i<; 2 -l)(3y 2 -l)log 



(to + l)(y - 1) 



(w-l)(y + l) 

2(y 6 w — w 6 y — 5y 4 w + 5w 4 y + 15y 2 w ~ I5yw 2 + 5w — 5y) 



6(y — w)(3wy + 1 



5A(y 2 - l) 2 (w 2 - 1) 



3 / 2 

+—w(w 
4 

2yw 5 



l)(3u> 2 - 2 

~l)(3y 2 -l)lo 
20yw 3 + 20y 2 w 2 



-wy(w 2 - 1) 



(y + !)(«;-!) 



_(l/-l)(io + l) 
20w 2 - 30yw - 15y 4 



5 + 3y 6 + 25y 2 



5(y 2 - l) 2 (w 2 - 1) 

^oo > 

^cs /-»cr 
^0 3 20 5 

o, 

3(y 6 - 5y 2 w 4 + 4yw 5 + 10y 2 w 2 - 5y 4 - 5w 4 + 10y 2 - 2§yw + 10w 2 ) 



5(y 2 - l) 2 (w 2 - l) 5 



0, 
-3 
KL4 



(y 2 + Ay + 5)(y - I) 2 



+ 2 (w - y) 



(y + 1) 2 

(Ayw 3 — 6y 2 w 2 



log 



V - 



1 



(y 2 - Ay + 5)(y + l) 2 



v u> — 1 J (y — l) 2 

lOu; 2 + y 5 w — 4y 2 w + 7yw + 6y 2 — Ay 3 



log 
-5-1 



(w 2 - l)(y 2 - l) 5 



(3.18) 
(3.19) 

(3.20) 

(3.21) 

(3.22) 
(3.23) 
(3.24) 
(3.25) 
(3.26) 

(3.27) 



We note that equations (3.22), (3.23) and (3.26) result directly from the initial conditions (3.12) and (3.13). They 
are consistent with equation (3.14). 
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FIG. 3. The source transfer functions T§ s (krj <C 1, dashed line) and (krj S> 1, solid line). 



Despite the complicated forms presented here, all these Green functions have simple asymptotic behaviors in the 
radiation- or matter-dominated regimes. Since we are more interested in the matter perturbations today and we know 



from equation (3.11) that Q c oc rj 2 oc a when r]/rj cq — ► oo, we can design a 'source transfer function' as 



T c (k;r))^ lim G c (k; r,, rj) 



Note that this is different from the definition of the standard CDM transfer function ( 2.66 ). Equations ( 3.1< 
then lead to the source transfer functions: 



rpc 



rpc 



T, 



10 



"■2oo 



rpc 



4.4 



(w 2 -l) (3w 2 -l)log 



w + 1 
w — 1 



6w 



TiCS 
n — 



1 lc 



2w 

5A(w 2 - 1) ' 

3, 9 . 9.n „ , f w — 1 
-( 3w 2 -2) + -w( W 2 -l)\og(— — 



rpCC 

2 io — 



rpCC 

1 2oo — 



5(w 2 - 1) 

rpCS 

4 



rpCl 

1 20 — ~ 



2 

5A 



w 



w 2 — 1 



- log 



IU + 1 

w — 1 



T cr — n 



5(w 2 - l) 2 
0, 



3 
10A 



2w 



log 



w + 1 
w — 1 



3.28) 



3.27) 



3.29) 

3.30) 

3.31) 

3.32) 
3.33) 
3.34) 

3.35) 
3.36) 

3.37) 
3.38) 



We plot these source transfer functions in Figures |^, |], and 0. They are now only functions of the initial time, but 
not of the final time. In the context of topological defects, the defect source was formed at rji <C ?7cq- Therefore 
it would be also interesting to investi gate th e asy mptotic behaviors of the source transfer functions with very early 
initial times. For r\\ <C ?7cq, equations (3.29)-(3.38) become: 
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^lco(i) — 2 ' ~~ ^2»(i) — T7^~~ 1°S i ^ioo(i) — r 2»(i) — > (3.39) 

3^-/r5cs _ mcc _ ~ cc _ a Vi mr _ ArpcT _ ^Veq fo An \ 

"y i 0(i) - J 10(i) - — 1 20(i) - A ~ J 10(i) ~ Ai 20(i) - 5^ > l d - 4U J 

where the subscript (i) denotes the condition 771 -C ?j e q- These asymptotic behaviors can be clearly seen in Figures 
|] and ||. We note that on sub-horizon scales, has a maximum at r\ ~ ?7 oq as seen in Figure S. Adding the fact that 
cosmic defects seed matter perturbations only on sub-horizon modes due to the compensation mechanism, it follows 
that the defect-induced matter perturbations are seeded mainly during the radiation-matter transition era. This is a 
generically different mechanism from inflationary models, in which matter perturbations are seeded during inflation 
in the deep radiation era when all the modes are well outside the horizon. Nevertheless, the defect and inflationary 
models both provide scale-invariant perturbations at horizon crossing, and these perturbations evolve similarly after 
horizon crossing. 



C. Degeneracy of the Green functions 



In principle we need ten Green functions (five for S c = 5l + <5f and five for S T = d] + Sf) in order to solve equations 
(3.1) and (3.2) by using the formalism ( |3.8| ) and (3.9). However, in addition to the dependence (3.14) by which we 
can reduce the effective number of the Green functions by two, there is another constraint we can invoke — the zero 
entropy fluctuation on super-horizon scales in the initial conditions, i.e. s = s = at 771 for modes k <C l/ffi (see 
eqs. pT8| and pTigQ ). Since the formation time 771 of the active source is normally so early that the condition k -C 1 /rji 
(and thus s = s = 0) is generally satisfied on the scales of our cosmological interest, we can rewrite equation (|3.8|) as 



5 n(Kv) = G3{k;f],i]i)S c (k,T]i) +Q± (fc; 77 , 77i)5 c (k,77i) , 



where 



■>N 



-,Nc 
H-7, 



-,Nt 
H-2 ) 



i = 3,4. 



From equations fl3.20| ), fl3.2l| ), (|3.24[ ) and j3.25[) , we can get 



fCC 

yic 



[5(t 



15u>V 



lby A 



-45w 2 y 2 + 2w 7 y - 10w 3 y - 6yw 5 - 50yw + 5 + 15u> 2 ) 



Using equation (3.14), we can also obtain Q\ = Q cs , so that 

r>c r>cs r*c /^cs 

These results yield the source transfer functions: 

1 



T: 



3oo 



rpCX 

1 lc 



T: 



3w 2 



30 



rpC rpCS 



5(w 2 - l) 2 

If the initial time is deep in the radiation era, i.e. rji <C 77 oq , we further have 



rpC rpCS 

1 40 — 1 



T 30(i) - 5A 2 V( 



cq -2 
— 2 K ??i . 



J 3oo(i) — J loo(i) 



40(i) 



0(i) 



2??eq -1 



^4oo(i) ~ J c 



• = log ipsa cx n 

0(1) 2»j cq & A?7i /: 



(3.41) 
(3.42) 

(3.43) 
(3.44) 

(3.45) 

(3.46) 
(3.47) 



where the last proportionality is only an approximation. Figure]^ shows the solutions of Tg and (= T-f^), while 
T| = Tq s and = are already shown in Figure |[ We note the the asymptotic behaviors indicated in equations 
( p6|) and ( |3~47l ) can be clearly seen in Figures || and @. Therefore, the original ten Green functions for solving S c 
and <5 r have now been reduced to four functions: two for S c (Q\ = C/ cs and and two for <5 r (Q\ = C? rs and £/ 3 ). 
We shall concentrate only on the solutions of <5 C , while leaving those of S r elsewhere To calculate <5f we need Q cs 
using equation (|3.9|); to calculate (5* we need = £ cs and £3 using equation ( [3.41 ). In solving Si, we note that £3 
transfers the initial perturbations of both matter and radiation [<5 c (?7i) + ^r(?7i)] to today, while G cs transfers the initial 
perturbations of their time derivatives [<5 c (?7i) + ^r(??i)] to the present. 
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FIG. 6. The source transfer functions T§ (dashed line) and T§ x (= Tf^; solid line). 



D. Solutions on intermediate scales 

With C/g and Q\ (— Q cs ) as the two basis Green functions, we can now work out the solutions on intermediate scales, 



using results derived in previous sections. In the matter era, Q\ <C Q1 on all scales so from equation ( 3.11 ) we know 
that Gio(k', VOiVi) — Qtooi^ Vo, Vi) when r/i ^> ?7 cq . This can be clearly seen from Figures^ and|(| In the radiation era, 

the perturbations [S c (r]i) + 8t(Vi)] or [5 c (r]i) + <5 r (?7i)] that were seeded well before the horizon crossing will evolve in 
the same way as in the standard CDM model due to the same zero entropy fluctuation initial condition. Therefore 
the solution interpolating between Ql {k\ 770, Vi) and G^^k; r)o, r/i) for r\\ -C ?7 cq will be the standard CDM transfer 
function. Thus we can write down a fit of the solution for the full gamut of k and r/i as 



g?(k; rjo, = Gingiva, m) + GioiVo^i) ~ &toc(Vo,Vi) T(k)I(k;rn), i = 3,4, 



(3.48) 



where 



T(k) 
I(k;m) 



(0.0534- 



2.75 ' 

l+3.83fc' 



ln(2e + O.llfe) 
1 + 30?7i 



30^(1 + ^)' 



(3.49) 
(3.50) 



and k is in units of 77^ (see equation ( All )). He re T (k) = T c (k, 770; ^bo = 0) is the standard CDM transfer function 
without baryons (modified from Ref. |l9|; see eq. | 2.66 | for t he de finition of T c (k, 770)), and I(k; 771) is a small correction 
near the horizon crossing to make the analytic solutions ( |3.48 ) fit the numerical results. For a given mode which 
is initially outside the horizon, the background contents of the universe compensate the defect source until horizon 
crossing. Therefore the detailed behavior of these Green functions near the horizon scale will affect the so-called 
compensation scale, beyond which no perturbations can grow. This means that the correction function I{k]rj\) in 
equation ( 3.48) ) actually plays an important role in getting the compensation scale right, and we shall discuss this 
IV C| . We have verified numerically for both and Q\ that the fit ( 3.4S ) is accurate within a 4% 



further in section 

error for any k and 771 (note that the initial conditions of Q\ and Q\ in the numerical verifications can be obtained from 
eqs. |3f] , [fU3| and fy§ ). Figure @ shows the numerical solutions of Q\ an d Q\ (=G CS ) wi thin a chose n do main of 
(k,r]i). It confirms the asymptotic behaviors indicated by equations in ( 3.45 ) (see also e qs. |fo.29 |, JT3C] and | 3.31 |), 
and plotted in Figures S and||. The asymptotic behaviors shown by equations ( 3.46 ) and ( |3.47 ) can be also marginally 
observed from Figure M. 
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k(h Mpc" 1 ) 

FIG. 7. The numerical solutions of G%{k\ rjo, r?i) (upper panel) and G%{k; r)o-,r)i) (=Q CB (k; rjo, r)\); lower panel). They both have 
been normalized to the scale factor today, ao/a oq . Each line has a different initial time 771, whose smallest and largest values 
are labeled in both plots. Successive lines have even logarithmic time intervals, and r/i is in units of rj cq . 
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k=k 




FIG. 8. Three domains on the (k, ?7i)-plane for the solutions of Green functions C/3 and Q\: Region I (k < fc cq = l/77 eq ), 
Region II (k > fc oq and k > l/r]i), and Region III (fc cq < k < l/r/i). These three regions are divided by the thick solid lines. 
Also shown are the rji = rjd (horizontal dashed line), and the k = k s (d) (vertical dashed line). 



Schematically, we can divide the (k, r7i)-plane into three regions for the solutions of Q\ (i = 3,4). As shown in 
Figure ||, these three domains are: Region I [k < k cq = l/rj e q), Region II [k > k cq and k > l/r/i), and Region III 
(fc eq < k < 1/rji). In Region I, the solution of G1 is Qf because the horizon crossing happens after ?7 q, after which 
Gfo = as argued before. In Region II, the solution is Gf^ because all modes in this region are inside the horizon 
all the time. We notice that Gi merges with Gf^ at the boundary of Regions I and II, where rji > r] oq . In Region 
III, the solution along the k direction is in the same form as the standard CDM transfer function. This is because 
modes with larger k cross the horizon earlier, so that their perturbations are suppressed after the horizon crossing 
for longer until ?7 cq . In addition, the solution along the 771 direction in Region III is in the same form as (? 4 C . This is 
because modes in this region are initially on super-horizon scales, and a given mode with different initial time 771 will 
experience the same amount of suppression resulting from the period between the horizon crossing and rj cq . Therefore, 
Regions I, II and III illustrate the intrinsic property of the solution ( 3.4S| ). 



E. The effect of baryons 



There is one important issue which we have not discussed — the effect of baryons. Prior to the photon-baryon 
decoupling at 77a, the CDM and baryons are dynamically independent. In this era, the photon-baryon fluid propagates 
as acoustic waves with a sound speed given by equation ( [2.11 ), preventing baryons from collapsing on small scales. 
Therefore there exists a sound horizon at the decoupling epoch d s (d) (hereafter simply the sound horizon) which is the 
distance such waves can travel prior to ?7d, and which is the largest scale at which the baryons can affect the evolution 
of density perturbations. It has been shown that inside the sound horizon d s (<j), not only are the CDM perturbations 
seeded before ?7d suppressed due to the presence of baryons (e.g. p4] , ^9|j2l| ) , b ut also the baryons themselves have an 
exponentially decaying power due to the Silk damping |22| (see also eq. [ 2.53 1), with acoustic oscillations due to the 
velocity overshoot p6Lp7|. After the decoupling, baryons evolve in the same way as the CDM does, so the matter 



perturb ations today can be obtained by linearly combining the CDM and baryonic fluctuations at rjd (see section II F 
and eq. | 2.65 |), and then evolving them to today. 

It follows that the baryonic effects tend to suppress the matter perturbations seeded before the decoupling epoch 
{v < f?d, see the horizontal dashed line in Figure^) and on scales inside the sound horizon (i.e. for k > fc s (d) ~ l/^ s (d)i 
see the vertical dashed line in Figure ||). The perturbations seeded after ij^ or on scales k < fc s ( d ) wi ll not be affected 
by the baryons. With this argument, we can impose a suppression factor on our current solution (3.48) to account 
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for the effect of baryons, i.e. the solution with the inclusion of baryons can be written as 

§f B \k; rj, 77;; ft, Q m0 , ^bo) = S-{k; r], m) B (k, r/i; ft, Cl mQ , ft B0 ), (3.51) 
where B(k, 77;; h, ^ m o, ileo) accounts for the baryonic suppression: 

[ 1, for 77i ^ r/ d or fc < fc s(d) ~ l/<4(d), 

where T(fc; ft, f2 m o, ^bo) is the usual standard CDM transfer function with the baryonic dependence. One accurate fit 
of T(k; ft, f2 m o, ^bo) is provided in Ref. [La]. We note that the ratio T(k; ft, f2 m o, ^bo) /T{k\ h, 1, 0) is unity outside the 
s ound horizon (k < fc s (d) = l/^s(d))i an d is less than unity inside the sound horizon. Referring to Figure |8j equation 



( 3.52 ) means that the value of B(k, rji] ft, Sl m o, Hbo) is less than unity in the region to the right and above the dashed 
lines, and is unity otherwise. We also note that in the low-f2 m o models, the sound horizon can be smaller than the 
radiation-matter equality horizon, i.e., it is possi ble th at fc s ( d ) = l/<4(d) > k cq f ^9| - In addition, there is a transition 



era (r/i < rjd) which is not included in equation (3.52). This is because in this era the baryonic effects do not fully 



operate as in the regime -C ?7d so that a good fit is not trivial to obtain. We have numerically verified equation 



( 3.52 ), although an accurate fit to the missing era 7fi < 7/d has yet to be found. 



F. Solutions in if / or A 7^ models 

The solutions we have obtained so far have assumed K = A = 0. For K =^ or A 7^ 0, the growing behavior of the 
CDM perturbations departs from that of a flat A = model only at very late times in the matter era (see later for a 
more detailed argument). This allows us to apply a universal suppression factor on to account for the effects of 

curvature or A: 

Sf ( fc 5 Vo, fj; h, Q m0 , Q B0 , tt A0 ) = Q mQ h 2 g(fl mQ , Q A0 )Sf {k; r] , fj; 1, 1, f2 B o, 0), (3.53) 
where k is in units of f2 m oft 2 Mpc -1 , and g(Q m o, Qao) is given by p8[ 

g(n m0 , n Aa ) = ^ T . (3.54) 

2 nzl - n A0 + (1 + n m0 /2)(i + n AO /70) 



In equation ( 3.53| ), the leading factor r2 m0 ft 2 results from the fact that the ratio of scale factors a /a cq is proportional 



to n m oh 2 and that the Green function Q^ B ^ = T^ B ^ao/a e q is proportional to this ratio. The factor g(Q m o,fl A o) 
accounts for the suppression of the linear growth of densi ty p erturbations in a K =/= or A-universe relative to an 
^mo = 1 and Q A p — universe |^8| (also verified in Ref. |p9|). The reason for k to have the unit il m oft 2 Mpc^Mn 
equation ( |3.53 ) is that the horizon size at radiation-matter equality r/ cq is proportional to (f2 m o/i 2 ) -1 (see eq. | All ] 



in Appendix |A|). 



For K ^ or A ^ 0, the extrapolation scheme (3.53) will be inaccurate when r\ is close to 770, i.e., when the 



background dynamics at 77 significantly departs from that of a flat A = model. Nevertheless, this extrapolation 
scheme is still appropriate for most models with active source for two reasons. First, in the context of cosmic defects, 
the power of matter perturbations on t he sc ales of our interest (k ~ 0.01-lftMpc" 1 ) is mainly seeded around 77 eq (see 
Figure [Io| and the discussion after eq. |3.40[ ). At this time, the curvature or A effects are negligible. Second, at late 
times when the curvature or A effects become important, these scales of our interest are already well inside the horizon 
so that any curvature terms in the perturbation equations can be neglected. Therefore, the only required change in 
the perturbation equations to account for the effects of curvature or A is simply to incorporate the correct background 
dynamics, and this involves only modifications in 0(77), ^0(77) and ^(77), whose solutions are given in Appendix [A] As 
can be seen in Figure the presence of curvature or a cosmological constant affects the background dynamics only at 
late times. More precisely, we verify that for (fl m o, ^ao) = (0-2, 0), (0.2, 0.8), (1, 0) and (2.0, 0), the largest observable 
scale for matter perturbations k w O.OlftMpc -1 corresponds to the horizon sizes at 77 ps 5, 5, 27, 54?7 cq respectively, 
whereas in these models the curvature or cosmological-constant domination occurs only at a much later epoch when 
77 > 770. At these moments (77 « 5, 5, 27, 54?7 cq ), the scale factor in the K ^ or A 7^ models departs from that in the 
fiat A = model only by less than one percent. Indeed, we have numerically verified that the extrapolation scheme 



(3.53) is accurate within a 5% error for 77; < 60/7 cq and fl A < 0.85 in A-models, for 77; < 20?7 eq and Q m Q > 0.2 in open 
A = models, and for 77; < 20077 eq and f2 m o < 2 in closed A = models. These ranges of cosmological parameters 
have apparently covered the values of our interest. 
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IV. IMPORTANT PROPERTIES 



With the Green-function solutions we have found, we can now analytically investigate some important aspects 
about the growth of cosmological matter perturbations. 



A. The standard CDM model 



First we investigate the relationship between our Green functions and the standard CD M tran sfer function, and 
thereby to justify the use of the standard CDM transfer function in the analytic solution ( 3.48 |). In the standard 



CDM model, there are no subsequent perturbations, so we have Sn = 6* N + 5% — S ! N . As discussed in equations ( 2.36| ) 
and (2.37), we also know that the CDM perturbations have a growing mode 5 c (k;rj) oc rj 2 on super-horizon scales 
(krj <C 1) for r\ 3> ffcq or rj <C rj cq . For the super-horizon modes in the radiation era and all modes in the matter era, 
this allows us to write 



5 c (k;r,)=A j (k)r, 2 , j = R,M, 



(4.1) 



where Aj is the coefficient of the growing mode in the radiation era (j =R: r\ <C ?7cq and kr\ <C 1) or in the matter era 
(j =M: r\ 3> ?7cq)- Thus using our Green-function solutions ( [3.41 ) and (3.48) with the initial conditions s = s = and 

5iy(k]r]i) = 2d n (k; rji) / rji as required by the adiabatic inflationary model, we can derive the standard CDM transfer 
function as 



A M 
An 



6 c (k; r)i)r) 2 
1 



eq 



4t7 c 2 



71 



4% 2 q 



A 2 r)?T< 



30(i) 



2A\X 



40(i) 



T(k) = -T(fc) 



(4.2) 



where we have used rj\ <C rj eq <C r\ and equations ( Al3j ), Q3.28 ), (3.46) and (3.47), and the last expression was obtained 



based on the formalism ( 3.4§| ). First, we note that the two terms involving T^i) and T£ ^ are equal, meaning that 



the two sets of initial perturbations [5 c (r){) + S r (r]i)] and [5 c (r)i) + <5 r (^i)] contribute equally to the present matter 
perturbations. Second, the T(k) here is nothing but the standard CDM transfer function which we have defined 
earlier. Third, the coefficient 2/5 in the final result of equation (4.2) is well known (e.g. |l7],Q), and here w e obt ained 
it using our Green-function solutions. This coefficient can be also obtained by first knowing from equation ( 2.45| ) that 
too is a constant on super-horizon scales (k <C l/?7), and then using its definition ( 2.38|) and equation (4.1) to compare 
its expressions for j =R, M. One will find Too = Aji/itG — 5Am/2ttG, which implies Am/Ar = 2/5 for k -C 1/rj. Thus 
the above derivation and result not only illustrate the relation between our Green functions and the the standard 
CDM transfer function T(fc), but also justify the use of T(k) in our formalism ( |3.48 ). 



B. Independence of the initial conditions 

One important problem for structure formation with causal seeds is to investigate how the source energy was 
compensated into the radiation and matter background when the seeds were formed at rj\. From the result ( 2.48| ) 



we know tha t the power spectrum of the pseudo energy r o should decay as fc on super-horizon modes. As argued 
in equation ( [2.49 ) , we can thus take tqo = as part of the initial conditions provided that the scales of interest are 



well outside the horizon initially. For similar reasons we can take s = s = 0, where s = 35 r /4 — <5 C . In addition, from 



equation (2.38) without baryons, we have 

too = e„ + JL (l) £ n N 6 N + ±^L. (4.3) 



8ttG \a ^ 4vrGa 

N- 



Since r o = is required at 771 , it follows that for a given O o(x, ryi), one can have different ways of compensating it 
into between (5/v and <5/y. It is thus vital to check the dependence of the resulting 8l(rj) on the way we compensate 
Ooo(x, rji) into the background initially. Consider the following two extreme cases, both satisfying Too = s = s = 0on 
super- horizon scales at 77;: 
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1. S c = 35 r / 4 = 0, 8 C = 3<5 r /4 = [— 47rG(a/d)0 o]i: Using equation ( 3.41 ), the normalized resulting initial pertur- 
bations can be calculated as 



Ai(r7,77i) 



Sc (k,rj) 
-47rG(a/a)e 00 ]i 



(4.4) 



2. 5 C = 38 T /4 = [-87rG(a/d) 2 e 00 /(4 - n c )] u S c = 3<5 r /4 = 0: Similarly we have 



A 2 (?7,?7i) 



{k,ri) 



. 2w(w 2 



-4^G(a/d)e 00 ]i 



'A(3w 2 + 1) 



(4.5) 



To see the difference in i5 c (k,rjo) today resulting from these two cases, one can calculate 



Ai 



2w(w 2 



1) Tl 



:Si) 



A(3w 2 + 1) T n c 



0. 



(4.6) 



where we have used equations ( 3.30 ), ( 3.45 ) and ( 3.48 ). This implies that no matter how the source 8oo(x, 771 ) is 
compensated into the background when it was formed (i.e. with any portions between 8n and 8n initially), it results 

in the same 8 C (fc, 770) today on scales which were outside the horizon at r]i. We note that this independence of the 
initial conditions was first numerically observed in Ref. flTfl , and here we have provided an analytic proof. 



C. Compensation and total matter perturbations 



With a complete set of Green functions for both initial and subsequent perturbations, we can now investigate 
the resulting total CDM perturbations and therefore the compensation mechanism in models with active source. 

Having seen the independence of the resulting S c (k, rj ) on the way the source energy is initially compensated into 



various background components, we can invoke equation ( |4.4| ) for $1, and equation (3.9) for 8^ to obtain <5 c (k, 770) = 
(5j(k, 770) + #c(k, 770). For a given mode at which krji <g; 1 initially, we have: 



5 c (k,T7o)=? c (k ) r7o)+^(k ) r / o) 

aim 



= 4vrG 



g cs (fc;77o,77 i )6oo(k,77 i ) 



g cs (k; 770, 77) e+(k, 77)^77 



>n 



8irGao 
5A 2 a C q 
8irGao 



-T{k)e Q0 (k,r)i) 



Vo a(fi) ~ 

T'(fc;77)^(e + (k,77)d77 



5A 2 



^{-T(fc)e o(k,77o)- 
a ca < 



</(> 



T'(k; 77)^|e+(k, 77) + r(fc)e 00 (k, 77) 



eh) 



where 



T'(k;fj) = 



Glo(k\<na,f\) 



(4.7) 
(4.8) 

(4.9) 

(4.10) 



and ^(fcjT/o,?)) is given by (3.48). The function T'(k;rj) is plotted in Figure ^. Here we notice that the quantities 
inside the outer most brackets in equations (4.8) and (4.9) are equivalent to nothing but the coefficient of the growing 
mode in CDM perturbations. Using equation (4J5), one can obtain the resulting perturbations 5 c (k, 770) by knowing the 
initial 0oo(k, 771) and integrating the evolution history of 6+ (k, 7 7). Hence this expression is convenient for numerical 
purposes. In addition, we see that the first term in equation (4.8) comes simply from the initial source energy, serving 
with an opposite sign to account for energy conservation. This is the so-called compensation. On the other hand, the 
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FIG. 9. The function T'(k;rj) (solid lines) and the standard CDM transfer function T(k) (the dashed line). Each solid line 
has different rj, whose highest and lowest values are labeled in units of 77 eq . Successive lines have even logarithmic time intervals. 



second term results from the subsequent evolution of 0+(k, 77), which actively creates the CDM density perturbations 
on sub-horizon scales (see later). This term also provides a way for defects to create non-Gaussianity. 

Alternatively, equation (4.9) provides a both physically and mathematically transparent way of interpreting how 
the perturbations are seeded by the source. First consider the integral term for a given mode k. When the mode is 
well outside the horizon, i.e. 77 <C 1/fc, T'(k; rj) equals T(k) by definition. Hence the two terms inside the inner brackets 
reduce to ©oi,i(k, fj)T(k) due to source stress-energy conservation (2.21). Since the power spectrum of 8oi,i(k, fj) falls 



off as fc 4 outside the horizon (see eq. |2.47|), we expect the quantity inside the brackets to be negligible until the given 
mode approaches horizon crossing. Near horizon crossing, Q 0ii (\&, fj) is no longer small, and T'(k;rf) starts departing 
from T(k) (i.e. T'(k;rf) w constant > T(k) oc fc~ 2 , see Figure g), so the two terms inside the inner brackets begin to 
contribute to the integral. This also explains why the correction function I(k; 77;) in equation fl3.48| ) is important in 
affecting the compensation scale. After horizon crossing, the significance of the two terms inside the inner brackets 
then depends on the subhorizon behaviours of their power spectra. 

As for the first term in equation fl4.9f ) , we see that for a superhorizon mode today, the integral in ([D]) is negligible 
as argued above so that only the first term contributes. It serves to give the opposite sign to the source energy so as 
to account for energy conservation on superhorizon scales today, and thus for the compensation at the present epoch 
770. On the other hand, if a given mode is well inside the horizon today, then the first term will be negligible provided 
that the source energy Ooo(k, 770) has a power-law fall-off inside the horizon, as it does for cosmic strings. Therefore 
in calculating CDM perturbations on scales of our interest, which are well inside the horizon today, the first term in 
equation ( |4.SD is negligible, so that it will not affect our compensation argument observed from the integral. 

This argument can be further strengthened by deriving the pseudo-energy today. From the definition of r o ( |4.3| ) 
and the final result of equation (|4.9|), one obtains 



?oo(k,77o) = (l-T(fc))G o(k,r7o) + 



T'(k; j))-^(e + (k, 77) + T(fc)6 o(k, rj) 



dfj . 



(4.11) 



From this result, one can clearly see that for super-horizon modes, T{k) is unity by definition so that only the integral 
survives. We have also seen from an earlier argument that on super-horizon scales, t he q uantity inside the square 
brackets is not hing but the 0oj,i(k, fj), which has a k 4 fall-off power spectrum (see eq. [ 2.47f ). It follows immediately 
from equation (4.11) that the pseudo-energy today, 100(770), nas a & 4 -decay power spectrum outside the horizon. This 
result confirms the super-horizon behavior of t 00 presented in equation ( 2.48| ). On the other hand, although (1 — T(k)) 
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is approximately unity for sub-horizon modes, the usual sub- horizon power-law decay in 0oo(k, t/q) (as in the case of 
cosmic strings) will still make the first term in equation (4.11) negligible inside the horizon. 

Thus we can see explicitly in a neat mathematical form how compensation acts on a given length-scale. From this 
analysis we can also see that the compensation scale is determined not only by the functions T'(k; rf) and T(k), but 
also by the properties of the source near the horizon scale. Once the detailed behavior of the source near the horizon 
scale is known, we can accurately locate the compensation scale using equation (4.9) or (4.11). We note that this 
result is different from the claim in Rcf. where multi- fluid compensation back-reaction effects were studied to 
show that the compensation scale arises naturally and uniquely from an algebraic identity in the perturbation analysis. 
Ref. |3l[] also investigated the compensation scale, and found constraints on the generation of super-causal-horizon 
energy perturbations from a smooth initial state, under a simple physical scheme. The compensation wavenumber was 
found to be constrained with fc c > 2i]~ 1 due to causality, depending on the behavior of the causal events. This result 
is not inconsistent with our finding above, where we further provide a quantitative way to locate the compensation 
scale for any given specific model. 



V. SUMMARY AND CONCLUSION 



In this paper we present a formalism which can be used to study the evolution of cosmological perturbations in the 
presence of causal seeds. In this formalism we invoked the fluid approximation in the synchronous gauge to model the 
contents of the universe, and assumed photon-baryon tight coupling until the last-scattering epoch to account for the 
baryonic effects. The approximation of instantaneous decoupling of photons and baryons was then employed at the 
last-scattering epoch. In particular, we demonstrated the accuracy of our formalism in the context of the standard 
CDM model, by comparing our results of density perturbations with those calculated from CMBFAST. 

We then derived the analytic solutions of matter density perturbations in a flat A = cosmology. The errors 
in Ref. | ft6| were corrected to yield a complete set of Green-function solutions for the super-horizon and sub-horizon 
modes (eqs. 3.3], [3.8], [3.9], 1 3 . 1 S ] — [ 3 .40 1 ) . The degeneracy among these Green functions was then found by comparing 
their initial conditions and employing the zero-entropy initial condition (eqs. [ 3.14|, 3.41]). This effectively reduces 
the number of the Green functions needed in the perturbation solutions (eqs. [3.3], [3.9|, [ |3.41 ]-[ 3.47| ]). With this great 
simplification, the solutions on intermediate scales were then easily found by the use of the standard CDM transfer 



function (eq. [3.48Q. This complete set of solutions were numerically verified to high accuracy. The baryonic effects 
were also considered (eq. J3.5l| ) . We then extrapolated these Green- function solutions to K ^ or A ^ models 
(eq. [ 3.53[ ]), with numerical justifications to high accuracy. 

Using these Green-function solutions, we investigated several important aspects of structure formation with causal 
source. We first demonstrated the relation between our Green functions and the standard CDM transfer function 



(eq. [4^]). Second we proved that the resulting matter perturbations to day is independent of the way the source was 
initially compensated into the background contents of the universe (eq. |4.Cl). With our Green- function solutions and 
the use of the pseudo-stress-ener gy t e nsor , we fi nally addressed the compensation mechanism in a mathematically 
and physically explicit way (eqs. [4 



1, | 4.11 |). In particular, the compensation scale was shown to be dependent 
not only on the dynamics of the universe, but also on the properties of the source near the horizon scale. Once given 
the detailed behavior of th e source near the horizon scale, the compensation scale can be accurately located using our 
Green functions (eq. [ 4.11 |). 

Although in the literature, there have been detailed treatments of theories with causal seeds, the formalism and its 
analytic solutions presented here will provide not only a physically transparent way for understanding the evolution 
of matter perturbations, but also a computationally economical scheme which is particularly pertinent when one 
needs to investigate the phase information of the resulting cosmological perturbations. Following the same line of 
development, we have been also working on the analytic solutions for radiation perturbations |2q| , which will be useful 
in computing the full-sky CMB anisotropies seeded by topological defects. Finally, we note that although we have 
been concentrating on investigating the perturbations with causal source, our Green-function solutions are completely 
general and therefore can be also applied to the study of models with acausal source. 
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APPENDIX A: COSMOLOGICAL BACKGROUND DYNAMICS 

With the discovery of the CMBR in 1964 |32|, the universe is believed to be mainly composed of not only matter 
but also radiation. After the discovery, several authors worked out the solutions in some FRW models with both 
radiation and matter 1 33 3f|. In this appendix, we aim to derive the general solution of FRW models, in the presence 



of both curvature and a cosmological constant. 

We assume that the universe is homogeneous and isotropic, and is filled with two fluids, radiation and dark matter, 
whose stress-energy tensors are also homogeneous and isotropic on average. We shall ignore the overall contribution 
of the stress energy from causal seeds like defect fields, because in general they are much smaller than the total 
energy density of radiation and matter. Thus in a FRW universe with only radiation and matter components that 
evolve independently and adiabatically, the scale factor a(rf) is determined by the unperturbed Einstein equation, or 
equivalently the Friedmann equation: 

•2 i is 2 87rG/3 m0 ao A 4 

a +Aa = (l + a) + — a, (Alj 

o o 

where a dot represents a derivative with respect to the conformal time r\, K is the curvature, p m is the matter energy 
density, A is the cosmological constant, and we have normalized a eq = 1. If we define 
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yixn 
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where H = a/a is the Hubble parameter, then we have from (Al) that Cl m + f2 r + fl A + FIk = 1 an d 
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^A0 __ 

Cl m0 8irGp m0 

We also notice that fi r o/^mo = ag <C 1. We define 
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-3K 



n 



nil) 



87rGp m0 ag 



2 ' 



2(V2- 1) 



£ = 



A'O 



c 



^mO a o 



(A2) 
(A3) 
(A4) 
(A5) 

(A6) 
(A7) 



where we note that B, C <C 1 due to ao > 1 and r2 m o ~ S^ao ~ ^ao according to the current observational results. 
Thus we can rewrite equation (Al) as 



^] = A 2 (l + a + Ba 2 + Ca 4 ), 

an; 



where 



A= — 



da 



VeqJo {l + a + Ba^ + Ca 4 ) 1 / 2 



A. 



(A8) 



(A9) 



Equation ( |A8| ) can then be numerically evaluated with certain choices of f2 m o, ^ao an d £^ao- Assuming three species 
of neutrinos and using p l0 = 2.0747 x lCP 51 GeV 4 J37j and the fact that at r) eq both the curvature and the cosmological 
constant terms are negligible in (Al), we obtain 



a = 23219 n m0 h 2 , 
Vcq = 16.310 (fi^/i^'Mpc, 
i cq = 3.4058 x 10 10 (a n0 fr 2 r 2 sec, 



where ry oq is in the units measured today. In certain cases, (A8) can be exactly solved: 

1. K = A = (i.e. n m0 = 1, n Aa = 0): 

a{rj) = A 2 i] 2 /A + Arj , 
t(rj) = A 2 if/12 + A v 2 /2, 

which give r/ cq = 3t cq / \/2. 

2. K < 0, A = (i.e. Q m0 < 1, ^ao = 0): 

a(r}) 



(A10) 
(All) 
(A12) 



(A13) 
(A14) 



2B 

«<*> = Jb 



cosh (AVBrj) + 2V B sinh {Aw Brf) - 1 
1 



cosh (Av Brf) 



sinh (AVBri) - — - 1 
2VB { " 2 



3. K > 0, A = (i.e. ft m o > 1, ft A0 = 0): 

1 



a(ri) = — cos (AV^Brj) - 2\^B sin (AV^Br]) - 1 
2B I 



1 



— 1 — -A77 

cos (AV^Bt?) + — == sin (AV^Bv) ~ -y - 1 



(A15) 
(A16) 

(A17) 
(A18) 



27 




We notice that at ea rl y times eq uatio ns (A15 - A16 ) and (A17-A18) reduce to equations (A13-A14). At late times 
equations (A13-A14), ( A15- Alq) and (A17-A18) give the asymptotic forms 



77 2 , K = A = 0, 

a(r)) cx { exp(AVBr]), K < 0, A = 0, 
1 - cos(Ay/^Br)), K > 0, A = 0, 



(A19) 



t 2 '\ K = A = 0, 

a(t) oc { t, K < 0, A = 0, 

1 - cos[2A(~B) 3 / 2 t], K > 0, A = 0. 



(A20) 



Figure |l0| shows some examples of these solutions. As we can see, the destinies of universes in different cosmologies 
diverge, although all have identical features around or before the radiation-matter equality t cq . This converging 
behavior at early times helps simplify the calculation of cosmological perturbations with causal source, since we know 
that this kind of perturbations are mainly contributed from the radiation-matter transition era. 
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